Transition frequency shifts with fine-structure constant variation for Yb II 
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In this paper we report calculations of the relativistic corrections to transition frequencies (q 
factors) of Yb II for the transitions from the odd-parity states to the metastable state 4/ 13 6s 2 2 Fy 2 - 
These transitions are of particular interest experimentally since they possess some of the largest q 
factors calculated to date and the 2 Fy 2 state can be prepared with high efficiency. This makes Yb II 
a very attractive candidate for the laboratory search for variation of the fine-structure constant a. 
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I. INTRODUCTION 

A discovery of acceleration of the universe (see, 
e.g., [H) is usually associated with the existence of the 
dark energy. The latter, according to the theories de- 
scribing cosmological evolution, may be a reason of vari- 
ations of the fundamental constants. There is an on-going 
discussion in the literature whether the fine-structure 
constant a could change during evolution of the universe 
or not. The Australian group reported in Q a nonzero 
result while other astrophysical groups do not confirm 
it [3, H[- But as it was argued in a more thorough 
analysis of the data used in [ij also leads to a nonzero 
result. New laboratory and astrophysical investigations 
are in progress. 

Laboratory studies of hypothetical variation of the 
fine-structure constant are based on the fact that tran- 
sition frequencies in atoms depend on aZ, where Z is 
the atomic number. Supposing that the modern value 
of a differs from its value in the earlier universe we can 
find relativistic transition frequencies shifts, determined 
by so-called q factors, according to 
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Most advantageous for these studies are the atoms and 
ions for which q factors of transitions between certain 
states significantly differ from each other. 

In @, it was proposed to use transitions whose q fac- 
tors are large and of opposite sign for laboratory measure- 
ments. In particular, it was shown that a good choice for 
an experiment would be Hg + because the q factor of the 
transition from the ground state to the low- lying 2 D 5 / 2 
state is very large and negative. There are two main 
considerations for the choice of a second transition for 
a comparison: 1) it should be convenient for an experi- 
ment; 2) its q factor should be small or positive because 
the measured quantity is proportional to the difference 
of q/u of the transitions being compared. 



At present the best laboratory constraint on the tem- 
poral variation of a of a/a = (—1.6 ± 2.3) x 10 _17 /year 
was obtained by Rosenband et al. in Ref. .7] by compar- 
ing the frequencies of the — ► 3 P§ transition in 27 A1 + 
and the 2 Si/ 2 — > 2 D$/2 transition in 199 Hg + . 

In this paper we propose to use Yb n for future exper- 
iments because there are several transitions whose wave- 
lengths can be synthesized with modern laser technolo- 
gies and whose q factors are very large and of opposite 
signs. 

Yb II has two low- lying metastable states 4f 14 5d 2 D 5 / 2 
(24333cm- 1 ) and 4/ 13 6s 2 2 F° /2 (21419 cm" 1 ) decay- 
ing to the ground state 4/ 14 6s 2 S , i/ 2 through electric 
quadrupole and electric octupole transitions, correspond- 
ingly. The q factors of the 2 D 5 / 2 and 2 F°, 2 states in re- 
spect to the ground state 2 S\j 2 were calculated in [f| H, [1] . 
For the former q w 10000 cm -1 and for the latter 
q w — 60000 cm -1 . In this work we carry out calcu- 
lation of the q factors for certain experimentally inter- 
esting odd-parity states with respect to the metastable 
state 4/ 13 6s 22 i 7 !^ /2 . Combining the results obtained in 
this work with the results presented in @, H, 0, we 
can easily find q factors ranging from — 60000 cm -1 to 
+75000 cm -1 associated with experimentally accessible 
transitions. This observation can potentially increase the 
sensitivity of the measurements for Yb n to a variation 
by more than a factor of two compared to all previously 
considered comparisons and makes this ion a very good 
candidate to establish a record constraint on a variaton 
in laboratory studies. 

The paper is organized as follows. Sec. [Til is devoted to 
the method of calculation of the properties of Yb II. In 
Sec. IIIII we discuss results and the experimental possibil- 
ities. Sec. IIVI contains concluding remarks. Atomic units 
(h = |e| = m e = 1) are used throughout the paper. 
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II. METHOD OF CALCULATION 

To find q factors we need to solve the atomic relativis- 
tic eigenvalue problem for different values of a or, respec- 
tively, for different values of a: from Eq. |T]). The value of 
x can be chosen somewhat arbitrarily, but two conditions 
should be satisfied. It should be sufficiently small to ne- 
glect nonlinear corrections and sufficiently large to make 
calculations numerically stable. Our experience shows 
that the choice of |x| = 1/8 allows us to meet both con- 
ditions. 

Thus, we have to calculate atomic frequencies ui± for 
two values x — ±1/8 of the parameter x. The corre- 
sponding q factor is given by 



q = 4(dj+ — LJ-). 



(2) 



TABLE I: The fifteen-electron CI calculations of the tran- 
sition frequencies to (in cm -1 ) for different basis sets. The 
Dirac- Coulomb Hamiltonian in the frozen-core approximation 
is used. The transition frequencies are calculated with respect 
to the metastable 4/ 13 6s 2 2 Fy 2 state. The experimental data 
from [13 ] are given in the last row. 



Config. 



Term J DHF [7sp5df] [8sp6df] Exp. 



4/ 13 ( 2 F°)6s 2 
4/ 13 ( 2 F 7 ° /2 )5d6s 



4/ 13 ( 2 F° /2 )5d6 S 
4f 13 ( 2 F° /2 )5d6s 



4/ 13 ( 2 F 5 ° /2 )5d6 S 
4/ 13 5d 2 

4/ 13 ( 2 F° /2 )5d6 S 
4/ 13 ( 2 F 5 ° /2 )5d6s 
4/ 13 ( 2 F 5 ° /2 )5d6s 
4/ 13 ( 2 F° /2 )5d 2 



The ground state configuration of Yb n is 
(Is 2 ... 4/ 14 6s). The configuration of the first excited 
state of Yb II is (Is 2 ... 4/ 13 6s 2 ). This is a metastable 
state because there is only weak (E3) transition from 
this state to the ground state. For experimental pur- 
poses this metastable state lives sufficiently long to be 
treated as a "ground" state. Since we are interested in 
calculation of q factors with respect to the metastable 
state with the open / shell, we have to treat Yb II as an 
ion with 15 electrons above closed shells. This makes 
calculations of Yb II rather complicated. 

In this paper we have carried out pure ab initio calcula- 
tions in the frame of the fifteen-electron CI method. The 
[ls 2 ...5p 6 ] electrons are treated as core electrons while 4/ ( F 7/2 )5rf 
4/, 6s, and 5d electrons are in the valence space. 

We started by solving the Dirac-Hartree-Fock (DHF) 
equations. The self-consistency procedure was done for 
(Is 2 ... 4/ 13 6s 2 ) configuration. After that, the 5d 3 / 2 an d 
5d 5 / 2 orbitals were constructed as follows: all electrons 
were frozen and one electron from the 6s shell was moved 
to the 5d shell, thus constructing the valence orbitals 
5d 3 /2 and 5d$/2 for the 4/ 13 5d6s configuration. 

In the next stage, we constructed virtual orbitals using 
the method described in [13, [H| and applied b y u s for 
calculating different properties of Fe I and Fe II [l2l [l3| . 
In this method, an upper component of virtual orbitals 
is formed from the previous orbital of the same symme- 
try by multiplication by some smooth function of radial 
variable r. The lower component is then formed using 
the kinetic balance condition. 

Our basis sets included s,p,d, and / orbitals with prin- 
ciple quantum number n < N\ t 2 which we designate as 
[NispN^df]. We carried out the calculations of energy 
levels, g, and q factors in a one-configurational (DHF) 
approximation and for two basis sets with (Aq = 7 and 
N2 = 5) and (N± = 8 and N2 = 6). Configuration space 
was formed with single and double excitations from the 
configurations 4/ 13 6s 2 , 4/ 13 5d6s, and 4/ 13 5d 2 . All re- 
sults which we discuss below are obtained in the pure 
Coulomb approximation (i.e, the Breit interaction was 
not included). 



3 [ll/2]° 
3 [ll/2]° 
3 [ll/2]° 
3 [5/2]° 
3 [7/2]° 
3 [9/2]° 

3 [9/2]° 
3 [7/2]° 
3 [9/2]° 
M7/2] 
1 [9/2]° 
3 [9/2]° 
3 [9/2]° 
3 [9/2]° 
3 [5/2]° 



1 [9/2]° 
3 [7/2]° 
3 [7/2]° 
3 [7/2]° 
3 [13/2]° 



7/2 














9/2 


5701 


8648 


7761 


8805 


11/2 


5810 


8835 


7985 


9144 


13/2 


6148 


9359 


8615 


10213 


7/2 


7493 


10564 


9831 


10561 


9/2 


8977 


12086 


11312 


11633 


7/2 


9603 


12578 


11737 


12076 


11/2 


15143 


13791 


12686 


13366 


9/2 


10613 


13713 


12940 


13600 


7/2 


10531 


13641 


12901 


13640 


11/2 


10834 


14645 


13575 


14413 


7/2 


18688 


17166 


15786 


16098 


7/2 


19610 


17782 


16384 


16923 


7/2 


15713 


18904 


18031 


18617 


9/2 


16599 


19734 


18866 


19499 


11/2 


17936 


21171 


20414 


21496 


7/2 


21092 


24366 


23634 


23019 


7/2 


22889 


25847 


23507 


23854 


9/2 


26117 


25510 


23633 


23916 


7/2 


20831 


23799 


22972 


24011 


9/2 


21626 


24719 


24412 


24751 


9/2 


23391 


26320 


24005 


24936 


11/2 


24585 


27315 


24679 


25129 



III. RESULTS AND DISCUSSION 

The results of the 15-electron CI calculations of the fre- 
quencies in respect to the metastable state 4/ 13 6s 2 2 Fy 2 
were obtained in the DHF approximation and for the ba- 
sis sets [7sp5df] and [8sp6df]. These results are given 
in Table [J 

As is seen from the table, the values found for the 
[7sp5df] basis set are in a good agreement with the exper- 
imental frequencies. For low-lying states this agreement 
is at the level of a few percent. For higher-lying levels 
the agreement is slightly worse but does not exceed 10%. 
The reason is that the configuration interaction is more 
significant for the high-lying states, as can be seen by 
comparison of the results obtained for the [7sp5df] and 
[8sp6df] basis sets. For the [8sp6df] basis set an agree- 
ment between theoretical and experimental values for the 
high-lying states belonging to the configuration 4/ 13 5d 2 
becomes much better. Note that the configuration space 
constructed for the [8sp6df] basis set included ~ 2 000 000 
determinants and the calculations of energies were rather 
lengthy and time-consuming for this basis so a further in- 
crease of the configuration space is impractical. 

At the same time, as it follows from Tables [TT] and 
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TABLE II: The fifteen-electron CI calculations of g factors 
for different basis sets. The Dirac-Coulomb Hamiltonian in 
the frozen-core approximation is used. The experimental data 
presented in the last row are taken from [l4| . 



Config. 



Term J DHF [7sp5df] [8sp6df] Exp. 



TABLE III: The fifteen-electron CI calculations of the q fac- 
tors (in cm -1 ) for different basis sets. The Dirac-Coulomb 
Hamiltonian in the frozen-core approximation is used. The 
q factors are calculated with respect to the metastable 
4/ 13 6s 2 2 F 7 ° /2 state. 



Config. 



Term J DHF [7sp5df] [8sp6df] 



4/ 13 ( 2 F°)6s 2 


2po 


7/2 


1.143 


1.443 


1.143 


1.145 




4/ 13 ( 2 F? /2 )5d6s 


3 [ll/2]° 


9/2 


0.934 


0.933 


0.933 


0.935 


4/ 13 ( 2 F 7 ° /2 )5d6s 


3 [ll/2]° 
3 [ll/2]° 
3 [5/2]° 


11/2 


1.122 


1.121 


1.121 


1.112 






13/2 


1.231 


1.231 


1.231 


1.230 






7/2 


1.377 


1.375 


1.371 


1.331 






3 [7/2]° 


9/2 


1.296 


1.286 


1.278 


1.264 






3 [9/2]° 


7/2 


1.027 


1.010 


0.995 


0.991 






1 [ll/2]° 
3 [9/2]° 
3 [7/2]° 


11/2 


1.093 


1.195 


1.131 


1.119 






9/2 


1.132 


1.140 


1.147 


1.158 






7/2 


1.048 


1.066 


1.084 


1.124 






3 [9/2]° 


11/2 


1.244 


1.143 


1.206 


1.214 






[7/2]° 


7/2 


1.130 


1.120 


1.121 


1.119 




4/ 13 ( 2 F° /2 )5d6s 


1 [9/2]° 


9/2 


1.079 


1.088 


1.089 


1.093 


4/ 13 ( 2 i ? 7°/ 2 )5d6 S 


4/ 13 ( 2 F 6 ° /2 )5d6 S 


3 [9/2]° 


7/2 


0.720 


0.720 


0.718 


0.720 


4/ 13 ( 2 F° /2 )5d6 S 


3 [9/2]° 


9/2 


0.984 


0.983 


0.982 


0.967 






3 [9/2]° 


11/2 


1.127 


1.128 


1.129 


1.115 




4/ 13 ( 2 F° /2 )5d6 S 


3 [5/2]° 


7/2 


1.002 


1.020 


1.045 


1.10 


4/ 13 ( 2 F° /2 )5d6s 


4/ 13 5d 2 


o 


7/2 


1.295 


1.297 


1.277 


1.18 


4/ 13 5d 2 


4/ 13 ( 2 F 5 ° /2 )5d6 S 
4/ 13 ( 2 F° /2 )5d6 S 


1 [9/2]° 


9/2 


0.945 


0.993 


1.000 


1.01 


4/ 13 ( 2 F° /2 )5d6 S 
4/ 13 ( 2 F 5 ° /2 )5d6 S 


3 [7/2]° 


7/2 


1.168 


1.158 


1.144 


1.150 


4/ 13 ( 2 F° /2 )5d6 S 
4/ 13 ( 2 F° /2 )5d 2 


3 [7/2]° 


9/2 


1.147 


1.092 


1.086 


1.10 


4/ 13 ( 2 F° /2 )5d6 S 

4f 13 ( 2 F° /2 )5d 2 

4/ 13 ( 2 F 7 ° /2 )5d 2 


3 [7/2]° 


9/2 


1.311 


1.310 


1.309 


1.29 


4/ 13 ( 2 F 7 ° /2 )5d 2 


3 [13/2]° 


11/2 


0.957 


0.955 


0.955 


0.97 



3 [ll/2]° 
3 [ll/2]° 
3 [ll/2]° 
3 [5/2]° 
3 [7/2]° 
3 [9/2]° 

3 [9/2]° 
3 [7/2]° 
3 [9/2]° 
1 [7/2]° 
M9/2] 
3 [9/2]° 
3 [9/2]° 
3 [9/2]° 
3 [5/2]° 



1 [9/2]° 
3 [7/2]° 
3 [7/2]° 
3 [7/2]° 
3 [13/2]° 



9/2 14522 

11/2 14843 

13/2 15396 

7/2 14563 

9/2 14736 

7/2 14876 

11/2 17113 

9/2 15372 

7/2 15518 

11/2 15270 

7/2 17892 

9/2 17243 

7/2 23242 

9/2 23795 

11/2 24610 

7/2 24243 

7/2 27078 

9/2 26185 

7/2 24611 

9/2 24432 

9/2 27489 

11/2 26804 



14458 
14813 
15548 
14481 
14690 
14761 
15910 
15462 
15568 
16988 
17698 
17541 
23388 
23968 
24951 
24568 
27201 
26562 
24180 
25180 
27711 
26816 



14118 
14496 
15371 
14167 
14402 
14469 
16573 
15306 
15396 
16015 
17691 
17577 
23072 
23650 
24767 
24439 
26910 
25918 
23976 
25521 
27638 
26607 



Mil where the results for g and q factors are presented, 
both these quantities are rather insensitive to the size of 
the configuration space. Using the experimental g fac- 
tors and comparing them with the calculated ones, we 
are able to identify properly the calculated energy lev- 
els. For the majority of the calculated levels there is a 
very satisfactory agreement between theoretical and ex- 
perimental g factors. This means that the configuration 
interaction is represented in a proper way for these states 
and also means that the q factors obtained for them are 
correct. 

Comparing the values of g and q factors obtained for 
the [7sp5df] and [8sp6df] basis sets, we see that almost 
all of them agree to each other at the level of few per 
cent. An exception is the levels with J = 11/2 belonging 
to the terms 1 [ll/2]° and 3 [9/2]°. These two states are 
closer in the [7sp5df] basis set approximation than in real 
experimental data so as to result in an artificial mixing 
of these levels and, consequently, to a change in their g 
and q factors. 

As is seen from Table IIIII the levels belonging to the 
4/ 13 5c£ 2 configuration have the largest q factors. This is 
not surprising because this configuration differs by two 
electrons from the configuration of the metastable state 
4/ 13 6s 2 . When the fine structure constant a tends to its 
nonrelativistic limit, one-electron energies of the 6s and 



5d electrons change in different directions, which leads to 
an increase of the q factors. Note in this respect that the 
high-lying levels which nominally belong to the 4/ 13 6s5<i 
configuration have (in reality) a rather large admixture 
of the 4/ 13 5<i 2 configuration. As a consequence their q 
factors are also large. 

We can estimate the accuracy of the calculated q fac- 
tors as a difference between the largest and smallest 
values (in each line) listed in Table IIIII for three basis 
sets. Such a conservative estimate shows that the ac- 
curacy of our calculations is not worse than 10%. In 
Table IIVI we present the recommended values of the q 
factors found in this work. For future reference, we also 
list in the table the q factors of the lowest-lying even- 
parity 4/ 14 5d 2 £> 3 /2,5/2 states and the 4/ 13 6s 2 2 F° /2 
state found in Refs. These values were obtained 

with respect to the ground state. Using the q factors 
found in this work with respect to the metastable state 
and q( 2 Fy 2 ) obtained in [1,0] with respect to the ground 
state, it easy to recalculate the q factors presented in Ta- 
ble [HI] with respect to the ground state. It can be done 
just by subtracting q( 2 Fy 2 ) from these q factors because 
for positive u> = u>i —lu2 the value of q is equal to (qi —52)- 

Experimentally one can search for a variation of a by 
comparing two frequencies of atomic transitions over a 
long period of time. Following the Ref. [8[ we can repre- 
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sent a measured quantity A(t) as 

w dt\LU 2 



0J\ UJ 2 
LU X UJ 2 ' 



(3) 



where u> = dio/dt. Taking into account Eq. ([T]) we can 
rewrite Eq. as follows 



A(t) 







- 




^2 / 


V "lab / 



(4) 



Of particular interest experimentally are narrow tran- 
sitions with large q values. Many of these ex- 
ist from the 2 Fy 2 metastable state to higher-lying 
states with wavelengths which can be synthesized 
with modern laser technologies. As is seen in 
Table IIV1 a very large A(t) can be expected if 
the frequencies of the 4/ 14 5d 2 D 3/2 (22961 cm" 1 ) - 
4 ji3 6s 2 2 F o (21419 cm" 1 ) and4/ 14 6s 2 S 1/2 (the ground 
state) — 4/ 13 6s 2 2 Fy 2 (21419 cm -1 ) transitions are com- 
pared. The q factor of the former (qi) is positive 
(qi « 71000 cm -1 ) and the q factor of the latter 
(q 2 ) is negative (q 2 « —60000 cm -1 ). The transition 
4/ 13 6s 2 2 Fy 2 ~4f 14: 5d 2 D b / 2 has nearly the same q value 
as the 4/ 13 6s 2 2 F y2 - 4/ 14 5d 2 D 3/2 transition and the 
necessary wavelength at 3.4 /im is somewhat easier to 
synthesize than the 6.5 /im to the 2 D 3 / 2 state. Another 
consideration is that the transition to 2 D 3 / 2 at 6.5 ^m 
is forbidden for even isotopes, but is weakly allowed 
through the hypcrfme interaction within odd isotopes. In 
either case, transitions from mp = — ► m' F = in odd 
isotopes are preferred due to the smaller (and quadratic) 
Zeeman shift to the transitions relative to the linear shift 
within the even isotopes. 

Substituting the q values into Eq. (|4]) for the two most 
sensitive transitions described above, and using the esti- 
mate |d/ai ab | < 10 -16 yr -1 @] we find 



J > 7/2 with energies less than 48000 cm -1 are most 
likely all metastable since there are no even-parity states 
with J > 5/2 below this energy. The g-factors are some- 
what lower, but this consideration may be outweighed 
by practical considerations such as available sources of 
narrowband laser light. The value of A{t) for a com- 
parison between 4/ 14 6s 2 Si/ 2 — 4/ 13 6s 2 2 Fy 2 and one 
of these transitions is still significantly greater than any 
other comparison considered to date. 

Preparation of the metastable state 4/ 13 6s 2 2 Fy 2 has 
already been performed in the laboratory by several 
groups including one of us (JRT). One of the simplest 
schemes to populate this state is from the 2 D 3 / 2 state 
at 22961cm -1 which is populated by spontaneous decay 
from 2 P°j 2 at 27062 cm -1 ; the upper state of the most ac- 
cessible laser cooling transition. Transitions from 2 D 3 / 2 
at 861 nm, 1062 nm and 1163 nm can all be used with 
varying efficiencies. For example, we have used a simple 
external cavity diode laser at 861 nm to prove the tech- 
nique. For any of these wavelengths, state preparation is 
limited by the 5 ms lifetime of the 2 D 5 / 2 state at 24333 
cm -1 . 



IV. CONCLUSION 

We have calculated relativistic frequency shifts (q fac- 
tors) for a number of transitions from excited states of 
Yb II to the metastable state 4/ 13 6s 2 2 F°, 2 state. We 



found that all these q factors are large (~ 10 4 cm 4 ) and 

1/2 



positive. The q factor for the transition 2 Si/ 2 — 2 Fy 2 is 



very large and negative. Because the 2 Fy 2 state is con- 
venient for an experiment and can be prepared with high 
efficiency, we conclude that Yb n is a very good candi- 
date for the laboratory search for possible variation of 
the fine-structure constant a. 



A(t) < 10 



-14 



yr 



(5) 



It is worth noting that a presence of transitions for which 
the q factors are very large and have opposite sign is 
favorable. First, it leads to increasing A(t) (as it follows 
from Eq. (TJJ) and, second, allows better control of some 
systematic errors that are not correlated with signs and 
magnitudes of the frequency shifts. 

It is also important to see that transitions from the 
4/ 13 6s 2 2 Fy 2 state to higher-lying states with J > 7/2 
possess more easily synthesized wavelengths. States with 
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TABLE IV: The recommended values of the q-factors (in cm found in this work with respect to the 4/ 13 6s 2 2 Fy 2 state and 
the q factors found in Refs. [(_,__,__] with respect to the ground state 4/ 14 6s 2 Si/ 2 . 



Config. 


Term 


J 


This work Ref. [9] Ref. [8J Ref. [6] 




2 D 


O/ — ■ 


inilS 19^89 




2 D 


5/2 


103Q7 11438 


Af 13 ( 2 F°\P.a 2 


2 p o 


7/9 


837^9 ^737 




3 [1 1 /2l° 

[11/ ZJ 


9 /2 


14100 


3 [ll/2]° 


11/9 


i /i ^nn 




3 [ll/2]° 


13/9 


i ^/inn 

lO^iUU 




3 [5/2]° 
3 [7/2]° 


7/2 


14200 




9/2 


14400 




3 9/2 ° 

L** / J 


7/2 


14500 




1 [11/21° 
3 [9/2l° 

L / J 


11/2 


16600 




9/2 


15300 




3 [7/2]° 


7/2 


15400 




3 [9/2]° 
M7/2] 
\V/2]° 
3 [9/2]° 
3 [9/2]° 


11/2 


16000 




7/2 


17700 
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